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Abstract 

We calculate the cross section and the transverse momentum distribution of inclusive 
J/ijj production in 77 collisions at CERN LEP2 within the colour singlet model and the 
/c-r-factorization approach, including both direct and resolved photon contributions. The 
unintegrated gluon distribution in the photon is determined, using the Blumlein's prescrip- 
tion for unintegrated gluon distribution in a proton. We compare our theoretical predictions 
with preliminary data taken by the DELPHI collaboration at LEP2. In addition, we present 
our predictions for the J/ip polarization properties. 



1 Introduction 

Heavy quark and quarkonium production at high energies has been vigorously studied 
in recent years from both theoretical and experimental viewpoint. At the modern collider 
conditions these processes are so called semihard ones [1-4]. In such processes by definition 
the hard scattering scale /1 ~ uiq is large compare to the Aq CD parameter but on the other 
hand /1 is much less than the total center-of-mass energy: Aqcd <C h <C \fs. The last 
condition implies that the processes occur in small x region: x ~ mq/^/s <C 1, and that the 
cross sections of heavy quark and quarkonium production are determined by the behavior 
of gluon distributions in the small x region. 
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It is also known that in the small x region it becomes necessary to take into account the 
dependences of the subprocess cross section and gluon distribution on a gluon transverse 
momentum fc T [1-4]. Therefore calculations of heavy quark and quarkonium production 
cross sections at HERA, LEP/LHC and other collider conditions are necessary to carry out 
in the so called A/p-factorization [2, 3] (or semihard [1, 4]) QCD approach, which is more 
preferable for the small x region than standard parton model. 

The /cT-factorization approach is based on Balitsky, Fadin, Kuraev, Lipatov 
(BFKL) [5] evolution equations. The resummation of the terms ln n (// 2 /AQ CD ), 

ln ra (/i 2 /AQ CD ) ln ra (l/a;) and ln n (l/x) in the /^-factorization approach leads to the 
unintegrated (dependent from q?) gluon distribution $(x, g^,/i 2 ) which determine the prob- 
ability to find a gluon carrying the longitudinal momentum fraction x and transverse momen- 
tum q T at probing scale fi 2 . In contrast with the usual parton model, to calculate the cross 
section of a physical process the unintegrated gluon distributions have to be convoluted with 
off mass shell matrix elements corresponding to the relevant partonic subprocesses [1-4]. 

Nowadays, the significance of the /^-factorization approach becomes more and more com- 
monly recognized. Its applications to a variety of photo-, lepto- and hadropro duct ion pro- 
cesses are widely discussed in the literature (see [6] and references cited therein). 

It is known that J/ip production at high energies is an intriguing subject in modern 
physics. It is traces back to the early 1990s, when the CDF data on the J/ip and T hadropro- 
duction cross section revealed a more than an order of magnitude discrepancy with theoreti- 
cal expectations. This fact has induced extensive theoretical investigations. In particular, it 
was required to introduce new additional production mechanism, the so-called colour octet 
(CO) model [7], where cc-pair is produced in the color octet state and transforms into final 
colour singlet (CS) state by help soft gluon radiation. Since then, the color octet model has 
been believed to give the most likely explanation of the quarkonium production phenomena, 
although there are also some indications that it is not working well. 

One of the problems is connected to the J / ifj photo- and leptoproduction data at HERA. 
For example, the contributions from the color octet mechanism to the J / ip photoproduction 
contradict the HI and ZEUS data for z distribution [8-1 1] . Also in leptoproduction case the 
shapes of the Q 2 , p\ and y* spectra are not reproduced by calculations [12] within the color 
octet model. The z distributions calculated later in [13] contradict the HI experimental 
data, too. 

Another difficulty of this model refers to the J / ip polarization properties in pp-collisions 
at Tevatron. If, as expected, the dominant contribution comes from the gluon fragmentation 
into an octet cc pair, the J/ip mesons must have strong transverse polarization [14-18]. 
This is in disagreement with the experimental data [19], which point to unpolarized or even 
longitudinally polarized J/ip mesons. 

Taking into account the above mentioned problems of color octet mechanism, we have 
considered the inelastic J ftp meson photo- and leptoproduction at HERA within the colour 
singlet model with /^-factorization [20]. Our theoretical results agree well with HI and ZEUS 
experimental data for all distributions without any additional J/ip production mechanisms 
(such as given by the CO model)(see also [21]). Also it was shown that the /c^-factorization 
approach gives a correct description of the J/ ip polarization properties in ep and pp interac- 
tions [20, 22] at the large transverse momenta px- 

Based on these theoretical results, in this paper we consider the inclusive J / ifj production 
at CERN LEP2 conditions within the colour singlet model with /c^-factorization. There 
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are several motivations for our study the J/ip meson production in 77 collisions (e + e~ — > 
e + e~ + J/ip + X). First of all, recently the DELPHI collaboration has presented preliminary 
data on the inclusive J/ip production in 77 collisions at CERN LEP2 [23], which wait to 
be confronted with different theoretical predictions. On the one hand it was shown recently 
that results obtained with account of the color octet contributions in the framework of 
the NRQCD do not contradict the DELPHI data [24]. On the other hand we know that a 
description of the inclusive J/ip production in pp collisions at Tevatron in the /^-factorization 
approach needs some significantly lower values for the color octet matrix elements to fit the 
Tevatron data [22, 25]. 

It is also known that the inclusive J/ip production in photon-photon collisions is dom- 
inated by the single-resolved process [26] and therefore it reveals the gluon structure of 
the photon. However the unintegrated gluon distributions in the photon $ 7 (x, q\, /i 2 ) are 
poorly known in contrast with the similar distributions in the proton and no attempts have 
been made to describe them until recently [27, 28]. Their knowledge is in particular neces- 
sary for the description of heavy quark and quarkonium production in 77 collisions within 
the semihard QCD approach. First application of the /c T -factorization approach for the 
case of resolved photons in heavy quark production is performed in [29] where the Kimber- 
Martin-Ryskin (KMR) [30] and Golec-Biernat-Wusthoff (GBW) [31, 32] prescriptions for 
the unintegrated gluon distribution in the photon were used. 

In the present paper we obtain in an independent way the unintegrated gluon distributions 
in the photon using the method proposed in Ref. [33] by J. Bliimlein for the gluon distribution 
in the proton. We calculate the cross section of inclusive J/ip production in 77 collisions 
including both direct and resolved photon contributions and compare them with preliminary 
experimental data [23] taken by the DELPHI collaboration at CERN LEP2. Additionally, 
we give our theoretical predictions for J/ip polarization properties at the LEP2 conditions. 

The outline of this paper is as follows. In Section 2 we present, in analytic form, the 
differential cross section for the inclusive J/ip production in the color singlet model with 
/cr-factorization and obtain the unintegrated gluon distributions $ 7 (x, gf., /i 2 ) in the photon. 
In Section 3 we present the numerical results of our calculations and compare them with the 
DELPHI data. Finally, in Section 4, we give some conclusions. 



2 Theoretical framework 

In this section we calculate total and differential cross section for inclusive J/ip produc- 
tion in 77 collisions within the color singlet model with /c^-factorization and obtain the 
unintegrated gluon distributions in the photon. 

2.1 Kinematics and cross sections 

In 77 collisions J/ip can be produced by one of the three mechanisms: a direct production 
(Fig. la), a photoproduction off a resolved photon (Fig. lb) and a by a double resolved 
process (Fig. lc). We will refer to the direct production, the once- and double resolved photon 
processes as direct, 1-res and 2-res ones respectively. The appropriate QCD motivated cross 
section the process e + e~ — > e + e~ J/ip X is given by 

da(e + e~ -> e + e~ J/ipX) = / f l/e (x 1 )dx 1 / f J / e (x 2 )dx 2 da('y-f -> J/ip X), (1) 
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Figure 1: Diagrams contributiong to inclusive J/ip production in 77 collisions. 

where we use the Weizacker- Williams approximation for the bremsstrahlung photon distri- 
bution from an electron [34]: 



with Q 2 min = m 2 e x 2 /(l - xf and Q 2 m£LX = (E6) 2 (l - x) + Q 2 min . Here x = E y /E e , E = E e = 
y/s/2, 9 = 32 mrad is the angular cut that ensures the photon is real, and y/s = 197 GeV [23]. 

For partonic cross section da(jj — > J/tpX) we will use expression obtained earlier [20] 
because the 1-res processes are analogous to those contributing to the lepto- or photopro- 
duction processes at HERA, and the 2-res ones are analogous to the J/ip hadropro duct ion 
subprocesses at Tevatron. In both cases the unintegrated gluon distributions in the proton 
should be replaced by ones in the photon. 

2.2 Unintegrated gluon distributions in the photon 

To obtain the unintegrated gluon distribution in the photon $ 7 (x, gf., /i 2 ) we apply the 
same method as for the gluon distibution in the proton according to the prescription given 
in [33]. The proposed method lies upon a straightforward perturbative solution of the BFKL 
equation where the collinear gluon density xG(x, /i 2 ) is used as the boundary condition. The 
unintegrated gluon distribution in a photon is calculated as a convolution of collinear gluon 
distribution iG 7 (i, /i 2 ) with universal weight factors: 

1 

$ 7 (x,g|,/i 2 ) = J p(r],q^,fi 2 ) ^-G 1 /i 2 ^ drj, (3) 



where 



a s 



as 



2 J ( 2 J a s ln(l/?7) ln(fi 2 /q 2 ) j , if q 2 < fj, 
I (2^asHlh)H4/^) , if Qt > V 



where Jo and Iq stand for Bessel function of real and imaginary arguments respectively, and 
as = 3a s /7T- In the capacity of the boundary condition we used expression for xG^x.n 2 ) 
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from the standard GRV set [35]. The parameter a$ is connected with the Pomeron tra- 
jectory intercept: A = 4asln2 = 0.53 in the LO and A = 4asln2 — Na 2 s in the NLO 
approximations, where N ~ 18 [36]. The latter value of A have dramatic consequences for 
high energy phenomenology. However, some resummation procedures proposed in the last 
years lead to positive value of A ~ 0.2 — 0.3 [37, 38]. Therefore in our calculations we will 
consider A as free parameter varying from 0.2 to 0.53 with a central value at A = 0.35. 




Figure 2: The unintegrated gluon distribution in the photon <& 7 (x, q\, M 2 ) as a function of longitudinal 
momentum fraction x (a) and transverse momenta squared q\ (b) of the gluon, for fixed values of x and q\. 
Curves 1 — 3 (Fig. 2a) correspond gluon virtuality q\ — 1, 10 and 50 GeV 2 as well as curves 1 — 4 (Fig. 2b) 
correspond longitudinal momentum fraction x = 10~ 4 , x = 10~ 3 , x = 10~ 2 and x = 10 _1 respectively. We 
set here = q\. 

We would like to note that intercept parameter A = 0.35 was obtained from the descrip- 
tion of pt spectrum of D* meson electroproduction at HERA [39]. This value of A also 
was used at the analysis of experimental data on inelastic J/ip photo- and leptoproduction 
at HERA [20], the bb quark production at Tevatron [40] and in the deep inelastic structure 
function F£, F£ and Fl description at small x region [41, 42]. 

In Fig. 2 the unintegrated gluon distribution in the photon (3) - - (4) (so called JB 
parametrization) is plotted as a function of longitudinal momentum fraction x (Fig. 2a) 
and transverse momenta squared q\ (Fig. 2b) of the gluon, for fixed values of x and q\. 
Curves 1 — 3 plotted on the Fig. 2a correspond gluon virtuality q\ = 1, 10 and 50 GeV 2 
as well as curves 1 — 4 plotted on the Fig. 2b correspond longitudinal momentum fraction 
x = 10 -4 , x = 10~ 3 , x = 10~ 2 and x = 10 -1 , respectively. 

Integrating expression (3) — (4) over the transverse momenta, we can obtain effective 
gluon distribution in photon and compare it to experimental data [42] taken by the HI 
collaboration at HERA. The results of such calculations are shown in Fig. 3. Curves 1 and 2 
correspond to the gluon density in photon from the standard GRV set [35] and the effective 
gluon density obtained from JB parametrization, respectively. One can see that both curves 
agree very well with the HI data and at x > 10~ 2 they practically coincide. 
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Figure 3: Effective gluon distribution in photon at Q 2 = 74GeV 2 . Curves 1 and 2 correspond to the gluon 
density from the standard GRV set [35] and the effective gluon density obtained from JB parametrization, 
respectively. Experimental data are from HI [43] collaboration •. 

3 Numerical results 

In this section we present the theoretical results in comparison with preliminary experi- 
mental data [23] taken by the DELPHI collaboration at CERN LEP2. 

There are three parameters which determine the common normalization factor of the cross 
section under consideration: J/ip meson wave function at the origin \I/(0), charmed quark 
mass m c and factorization scale \i. The value of the J/ip meson wave function at the origin 
may be calculated in a potential model or obtained from the well known experimental decay 
width r(J/V> -> In our calculation we used |^(0)| 2 = 0.0876 GeV 3 as in Refs. [15, 

20]. 

Concerning a charmed quark mass, the situation is not clear: on the one hand, in the 
nonrelativistic approximation one has m c = m^/2 = 1.55 GeV, but on the other hand 
there are many examples when smaller value of a charm mass is used in J/ip production 
processes at high energies, for example, m c = 1.4 GeV [12-15, 24]. However, in our previous 
papers [20] we analyzed in detail the influence of the charm quark mass on the theoretical 
results. We found that the main effect of a change of the charm quark mass connects with 
the final phase space of the J/ip meson, and in the subprocess matrix elements this effect is 
negligible. Taking into account that the value of m c = 1.4 GeV corresponds to the unphysical 
phase space of the J/ip state, in the present paper we will use the value of the charm mass 
m c = 1.55 GeV only. 

Also the most significant theoretical uncertanties come from the choice of the factorization 
scale up and renormalization one /xr. One of them is related to the evolution of the gluon 
distributions $ 7 (:r, q?, Hp), the other is responsible for strong coupling constant as(^%)- As 
often one done in literature, we set them equal /ip = //r = /i. In the present paper we used 
the following choice /i 2 = q\ as in Refs. [4, 20]. 

The results of our calculations are shown in Fig. 4 — 6. Fig. 4 displays the J/ip meson 
transverse momentum distribution in comparison with DELPHI preliminary data [22] at 
^/s = 197 GeV and \y^\ < 2. The solid lines 1 and 2 correspond to the standard parton 
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Figure 4: The J/if> meson transverse momentum distribution in comparison with DELPHI preliminary 
data [23] at y/s — 197 GeV and < 2. The solid lines 1 and 2 correspond to the standard parton model 
(SPM) calculations at the leading order approximation with the GRV (LO) gluon density (at m c = 1.55 GeV) 
and the fey-factorization results with the JB unintegrated gluon distribution (at m c — 1.55 GeV and A = 
0.35). The shaded band of the fey-factorization predictions connected with the value of pomeron intercept 
A: the lower border corresponds to A = 0.2 and the upper one corresponds to A = 0.53. 

model (SPM) calculations at the leading order approximation with the GRV (LO) gluon 
density (at m c = 1.55 GeV) and the /cy-factorization results with the JB unintegrated gluon 
distribution in photon (at m c = 1.55 GeV and A = 0.35). The shaded band is the fcy- 
factorization predictions connected with the value of pomeron intercept A: the lower border 
corresponds to A = 0.2 and the upper one corresponds to A = 0.53. One can see that the 
behavior of the transverse momentum distribution within the experimental and theoretical 
uncertanties can be explained by the color singlet model with fey-factorization 3 , while the 
collinear approximation prediction are lower than data by an order of magnitude. 

We would like to note a some difference in the shapes between curves obtained using 
the fey-factorization approach and the colliner parton model. This difference manifests the 
Pt broadening effect which is connected to the initial gluon transverse momentum and it is 
usual for the fey-factorization approach [20, 39, 40]. 

Fig. 5 show the transverse momenta and rapidity distributions of the inclusive J/ if) meson 
production at yj~s = 197 GeV, where all contribution (direct production, 1-res and 2-res 
processes) are plotted separately. Curves 1 and 2 and shaded bands are the same as in Fig. 4. 
One can see that 1-res contribution dominates, while the direct and 2-res contributions are 
only a small corrections to the 1-res one. This is contrast to the inelastic J / if) production case 
at HERA, where the resolved photon contribution is an small fraction of the direct one. Also 
one can see that 2-res contributions is more sensitive to various theoretical uncertainties. 

As it mentioned above, the main difference between fey-factorization and other approaches 
connects with polarization properties of the final particles. The " nonsense" polarization of 

3 In a Ref. [24] it was claimed about some evidence for color octet mechanism of NRQCD from these DELPHI experimental 
data [23]. However in a recent paper [44] the Belle collaboration has published results on prompt J/tjj and double cc quark 
production in e + e _ -annihilation, where the additional cc pair fragments into either charmonium or charmed mesons. It was 
shown that J/ip signal predicted in color octet model is not observed and that large fraction of prompt J /if) events is due to 
the e + e — — * J/ipcc process [44], We would like to note it is possible that in the DELPHI data at LEP2 the contributions of 
J/ipD*^ associated production are presented too. The calculations of these processes are in progress. 




Figure 5: The transverse momenta and rapidity distributions of the inclusive J ftp meson production at 
y/s = 197 GeV, where all contribution (direct production, 1-res and 2-res processes) are plotted separately. 
Curves 1 and 2 are the same as in Fig. 4. 

the initial BFKL gluons should result in observable spin effects of final J/ip mesons [45]. In 
the present paper we calculate the p\ dependence of the spin aligment parameter a: 

2 d<jjdp\ - 3 do L jdp\ , , 

a[PT >~ da/dp* +da L /dp* ' 1 j 

where is the production cross section for the longitudinally polarized J/ip mesons. The 
parameter a controls the angular distribution for leptons in the decay J ftp — » fi + yT (in the 
J/i/j meson rest frame): 

dT{J/ip -> ^+ //-) 2 

~ 1 + a cos 0. (6) 

d cos 

The cases a = 1 and a = — 1 correspond to transverse and longitudinal polarizations of the 
J/ip meson, respectively. 
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Figure 6: The y\ dependence of the spin aligment parameter a for the inclusive J/ip meson production at 
LEP2. Curves 1 and 2 are the same as in Fig. 4. 

The results of our calculations are shown in Fig. 6. Curves 1 and 2 are the same as in 
Fig. 4. As in inelastic J /if) photoproduction at HERA [20], we have large difference between 
predictions of the leading order of collinear parton model and the /c^-factorization approach. 
It was shown that account of the color octet contributions do not change a predictions of kr- 
factorization approach for the J /if) polarization properties in pp and ep collisions (see [22]). 
Therefore experimental study the polarized J /if) production at CERN LEP2 will be an 
additional test of BFKL gluon dynamics. 



4 Conclusions 

In this paper we considered the inclusive J / if) meson production at CERN LEP2 within the 
colour singlet model and the /c^-factorization approach, including both direct and resolved 
photon contributions. The unintegrated gluon distribution in the photon is determined, using 
the Blumlein's prescription for unintegrated gluon distribution in a proton. We compared the 
theoretical results with preliminary experimental data taken by the DELPHI collaboration at 
CERN LEP2. We find that the fc^-factorization results (in contrast with the SPM ones) with 
the JB unintegrated gluon distribution in the photon agree with data within the experimental 
and theoretical uncertanties. Finally, it is shown that experimental study of a polarization 
of J I if) meson at LEP2 should be additional test of BFKL gluon dynamics. 
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